Biophysical Reports

The mechanical properties of red blood cell (RBC) membranes are critical to their function in oxygen delivery, and changes to these properties as RBCs age affect their journey around the circulatory system, including clearance by the spleen. Such changes can also have significant health effects,including in cardiovascular disease and on the interactions between RBCs and malaria parasites. The function of RBCs requires them to be very soft, which together with the small size of the cells brings the energy required for measurable deformation of the membrane into the range of typical thermal energies. Red Blood Cells can be observed to flicker under optical microscopy, and these shape fluctuations can be quantified and used to obtain key biophysical parameters such as the tension and bending modulus of the membrane. Typically the shape of the cell's equator is extracted, and the mean power spectrum is obtained by time-averaging the power present in the normal modes of the thermal fluctuations. Flickering spectroscopy has been used extensively on RBCs, but has so far been very low throughput and with some technical limitations. Here we address issues related to active versus passive fluctuations, focusing and optics, camera exposure and sampling, and fast contour detection. In combination with an automated imaging system it is possible to measure thousands of cells in one day, with no user input. We validate this new pipeline by chemical modification of RBC mechanics, and by comparison with simulated fluctuations including each confounding effect. The methods and codes for this robust flickering analysis will allow consistent measurements across labs.

Fluorescence lifetime imaging microscopy provides quantitative, concentration-independent contrast for probing molecular interactions, biochemical environments and cellular physiology. However, the requirement to acquire sufficient time-resolved photon statistics makes FLIM inherently slow, limiting its application to dynamic biological processes. In live-cell applications, including calcium signalling and vesicular trafficking, acquisition times can exceed the timescale of the underlying biology, causing temporal averaging, motion artefacts and loss of transient information. Methods that increase FLIM acquisition speed while preserving quantitative lifetime accuracy are therefore required. We report on a high-speed, compressive, multiphoton fluorescence lifetime imaging technique (TomoFLIM). Two-photon fluorescence is excited using a line focus projected tomographically across the sample, while time-tagged fluorescence is acquired using time-correlated single-photon counting. Time-resolved fluorescence data are reconstructed using Lucy-Richardson deconvolution followed by a computationally efficient centre-of-mass method lifetime estimator. In addition, TomoFLIM Net, a physics-informed neural-network, directly reconstructs fluorescence intensity and lifetime from compressed time-resolved tomographic data. TomoFLIM was benchmarked against raster-scanned fluorescence lifetime measurements using calibrated fluorescence lifetime beads and biological specimens. We demonstrate imaging at compression ratios exceeding 90%, with Pearson correlation coefficients above 80% relative to reference images. A raster-scanned FLIM dataset acquired in 60 s was reproduced using TomoFLIM in 3.75 s, representing a 16-fold increase in frame rate and equivalent reduction in accumulated dark counts. TomoFLIM Net recovered distinct experimental bead lifetime populations, demonstrating a direct route from compressed measurements to quantitative lifetime maps. TomoFLIM therefore offers significant potential for rapid live-cell imaging of dynamic biological processes, including deep within turbid biological specimens.

Diffusion within the cytoplasm is fundamental to numerous biological processes. Fluorescence recovery after photobleaching (FRAP) is one of the most common method for quantifying molecular diffusivity in living cells using standard laser scanning confocal microscopy (LSCM). However, accurately measuring fast cytoplasmic diffusion (typically >10 m^2/s) is challenging due to rapid recovery kinetics, weak signal-to-noise ratios, post-bleach signal artifacts, and spatial restrictions affecting normalization. While individual challenges have been addressed in specific contexts, a simple and robust framework to quantify cytoplasmic diffusivity remains elusive. Here, we present a FRAP methodology specifically designed to overcome these obstacles. By utilizing the Gaussian function -- the impulse response (ImpRes) of the diffusion equation in an infinite medium -- our approach leverages the full spatiotemporal dataset through a single-equation three-parameter fitting procedure, thus releasing restrictions to small regions of interest and arbitrary initial time-points. The methodology was validated on three datasets of increasing complexity: in silico simulated recovery profiles, in vitro data from FITC-dextran in glycerol solution, and live-cell imaging of free cytoplasmic GFP. Systematic comparison with existing models demonstrates that the ImpRes approach significantly reduces sensitivity to noise and imperfect fluorescence normalization, while remaining robust against short-term biases, such as transient probe photo-activation. Given its robustness under realistic experimental conditions and its ease of implementation, the proposed FRAP methodology provides a reliable tool for quantitative cytoplasmic analysis.

Finite-bandgap charge-based FET biosensors are known to be insensitive to electronic tunneling noise in the semiconductor and interference from weakly-charged proteins in physiological samples. However, their signal is still corrupted by pH fluctuations and condensing counterions in the bulk electrolyte. Here, we demonstrate that electrolyte noise is screened by the Stern layer beyond the linear Debye-Huckel limit, when the interfacial potential drop is larger than the thermal voltage. Using a fully nonlinear Gouy-Chapman-Stern (GCS) theory, a universal calibration curve is derived and shown to collapse the low-noise FET signals for miRNA, extracellular vesicles (EV), and protein, captured either directly or through silica and magnetic nanoparticle charge reporters, over eleven decades of analyte concentration. The theory involves the dissociation constant KD and a gain parameter {chi}, which represents the competing effects of the reporter charge and buffer Debye screening that must be large to produce the low-noise GCS condition. The limit of detection (LOD) is shown to be KD/{chi}, which can be optimized by varying the probe density, charge reporter, and buffer solution. An optimized hybridization buffer then yields {chi} ~ 22 (LOD ~1 pM) for miRNA detection, and multivalent capture reduces KD by four orders of magnitude for EV detection (relative to protein detection) with {chi} ~ 34 (LOD ~10 fM).

Dynamic conformational changes in intracellular domains fundamentally affect the function and pharmacology of many membrane proteins. For example, sodium influx through the cardiac voltage-gated sodium channel (NaV1.5) is rapidly terminated through conformational changes that result in pore closure, a transition known as inactivation. Inactivation involves NaV1.5 intracellular regions, particularly the DIII-DIV linker containing the IFM particle (Isoleucine-Phenylalanine-Methionine) and its dysfunction is a major cause of cardiac arrythmias. However, the conformational changes involved in inactivation and their modulation by auxiliary proteins and clinically used drugs remain incompletely characterized, partly because live-cell, site-specific labeling of intracellular regions with small fluorescent dyes remains challenging. Here, we combine live-cell protein semi-synthesis with voltage-clamp fluorometry (VCF), to track intracellular conformational dynamics of the cardiac sodium channel NaV1.5 and monitor their voltage dependence and kinetics in real time. We identify intracellular conformational changes involved in both fast and steady-state inactivation of NaV1.5 and show that both lidocaine and auxiliary proteins affect the kinetics of conformational changes of the DIII-DIV linker. Our work establishes the combination of protein semi-synthesis and voltage-clamp fluorometry as a powerful approach to dissect intracellular conformational changes in membrane proteins.

Advances in single-photon sensitive detectors are rapidly expanding the adoption of photon-counting fluorescence microscopy. Under the Poisson photon-counting statistics, iterative Richardson-Lucy (RL) -type algorithms are statistically optimal for image deconvolution but suffers from a fundamental semi-convergent behaviour: prolonged iterations inevitably amplify noise, requiring heuristic early stopping or regularisation typically based on assumptions about object morphology. Here we present a morphology-agnostic regularisation framework for RL deconvolution that exploits independent noise realisations instead of structural object priors. Preserving the Poisson statistics of the acquisition process, we formulate the Regularized by Noise (RbN): a regularized variational framework and its corresponding iterative minimization algorithm that exploit the statistical consistency of independent noise realisations to distinguish reproducible image features from stochastic noise. The regularisation strength is selected automatically using the Poisson residual whiteness principle, resulting in a fully data-driven reconstruction without heuristic parameter tuning. Photon-timing-resolved systems naturally provide the independent noise realisations exploited by our framework, whereas computational photon splitting provides a statistically equivalent implementation for photon-counting systems. We validate the approach experimentally using photon-timing-resolved confocal microscopy and image scanning microscopy across eight morphologically distinct subcellular targets, and further demonstrate its applicability to photon-counting microscopes through computational photon splitting. Across imaging modalities, detector technologies, biological structures and signal-to-noise regimes, our method eliminates RL semi-convergence, removes sensitivity to the stopping criterion, and consistently outperforms conventional RL while preserving fine structural detail. More broadly, our results establish independent noise realisations as a general source of morphology-agnostic regularisation. Although demonstrated here for SPAD-based laser-scanning microscopy and Poisson statistics, the underlying principle could be extended to other imaging modalities and, more generally, to statistical inverse problems through appropriate noise-specific formulations.

Weak, reversible interactions underpin biomolecular recognition, and single-molecule fluorescence (smF) imaging techniques can provide unprecedented insight into those biological processes. Unfortunately, such studies often require micromolar concentrations of fluorophore-labeled biomolecules, which is beyond the accessible range of conventional smF microscopies. Here, we describe a surface-functionalization method based on cloud-point polyethylene glycol (PEG) grafting that enables widefield smF microscopy measurements at micromolar concentrations without the use of nanophotonic devices. Using conventional total internal reflection fluorescence (TIRF) microscopy, we detected single-molecule fluorescence resonance energy transfer (smFRET) from surface-tethered, donor-labeled target molecules with up to 8 micromolar concentrations of freely diffusing, acceptor-labeled analyte molecules in the background--two orders of magnitude higher than typical studies in the literature. Weak, DNA-hybridization and protein-RNA binding equilibria were measured across micromolar range titrations. Together with advances in high-background data analysis, the robust method presented here enables kinetic and thermodynamic analyses of weak biomolecular interactions, especially those limited by nonspecific adsorption and high fluorescence backgrounds, using only standard smF instrumentation.

Quantitative fluorescence imaging techniques such as fluorescence lifetime imaging microscopy and hyperspectral imaging infer molecular contrast from photons distributed across spatial pixels and temporal or spectral channels. In the few-photon regime, however, conventional pixel-wise analysis discards the spatial relationships imposed across neighboring pixels by the microscope point-spread function (PSF). Here we show that this spatially distributed information can be recovered without prior knowledge of emitter positions, spatial support or component assignments. We introduce SPOOL (Spatially Pooled Optical Observation Likelihood), a training-free Poisson inverse framework that jointly recovers source-space amplitudes and quantitative contrast by combining the PSF with temporal-decay or spectral-response dictionaries. For an isolated source, the attainable precision gain is governed by a dimensionless optical quantity: the PSF width expressed in detector pixels. The predicted gain therefore scales with optical sampling rather than with the physical origin of the contrast. The model predicts that lifetime-precision gain scales approximately linearly with the number of pixels spanning the PSF full width at half maximum, a scaling reproduced by Monte Carlo simulations. At one detected photon per foreground pixel, the reconstruction reduces lifetime dispersion sixfold in fluorescent-bead experiments and decreases the lifetime root-mean-square error relative to a high-photon reference from 1.19 to 0.45 ns in dual-labeled cells. The same framework transfers unchanged to hyperspectral imaging, recovering spectral contrast from generic emission bands without prior fluorophore spectra.

Nanopore sensors have transformed single-molecule analysis, enabling real-time detection of biomolecules with unprecedented resolution. Understanding and modelling noise in nanopore sensing is essential to unlocking their full analytical potential. Although mechanistic models of nanopore noise exist, a rigorous statistical framework that characterises noise structure across diverse experimental conditions is still lacking. Here, we present a comprehensive empirical analysis of nanopore noise using twelve diverse experimental conditions acquired with different nanopores, analytes, and electronic recording systems. By employing multiple, mutually reinforcing statistical methods, we establish that nanopore noise is well described by a composite fractional model comprising multiple fractional Gaussian noise and fractional Brownian motion components. This finding is supported by confirmation of Gaussianity, rigorous interpretation of second-order exponents, multifractal analysis, and evidence that transient deviations are attributable to deterministic or experimental artefacts rather than alternative stochastic mechanisms. We further resolve quasi-deterministic structures, such as baseline trends and structural breaks, and demonstrate their separability from the underlying stochastic profile. Collectively, these results provide the first formal statistical validation of a composite fractional model for nanopore noise and establish a principled basis for noise reconstruction, realistic simulation, and statistically grounded algorithm design for nanopore sensing applications. This study provides the foundation for the creation of high-fidelity performance evaluation datasets for signal processing, and the controlled generation and augmentation of data for machine learning while maintaining statistical fidelity to experimental conditions.

With its unique spatiotemporal resolution at the single-molecule level, Forster resonance energy transfer (FRET) is a powerful tool for investigating biomolecular conformational dynamics and function. However, early photobleaching remains a key limitation restricting the achievable observation time and therefore the information gain and application range of single-molecule FRET. Here, we establish nanophotonic DyeCycling, which overcomes photobleaching by reversible fluorophore binding and efficient background suppression using zero-mode waveguides. Nanoscale conformational changes are sensitively resolved from milliseconds up to the hour range, enabling robust kinetic analyses at the single-molecule level. As a result, outlier molecules reflecting static heterogeneity and time-dependent kinetics within individual molecules reflecting dynamic heterogeneity are reliably resolved. The versatility of DyeCycling is demonstrated using DNA and protein systems. Altogether, nanophotonic DyeCycling extends smFRET observations beyond the conventional photobleaching limit, revealing previously inaccessible kinetic effects and providing deeper insight into biomolecular systems with multiple states and rates.

Fluorescent proteins are widely used as quantitative reporters of protein abundance in living cells. Here we show that this relationship can break down under translation inhibition: the cellular ribosomal abundance can increase while fluorescence intensity remains unchanged or even decreases. Using fluorescent reporters of ribosome abundance in Escherichia coli and Bacillus subtilis, we find that cellular fluorescence intensity quantitatively tracks ribosome abundance, inferred by RNA-to-protein ratio, under nutrient-limited growth but that is no longer true during chloramphenicol treatment. In E. coli, the discrepancy occurs with transcriptional and translational reporters and with both ribosomal protein and ribosomal RNA promoters. To distinguish changes in fluorescence output per reporter from changes in reporter abundance, we constructed a fluorescent protein-LacZ dual reporter in which fluorescence and an independent enzymatic estimate of reporter abundance are obtained from the same protein. Under translation inhibition, LacZ-derived reporter abundance increases whereas fluorescence intensity remains approximately constant, showing that fluorescence output per unit reporter abundance decreases. In a companion study, Bakshi and colleagues demonstrate that the effect extends to other fluorescent proteins and other translation inhibitors, indicating that it is not specific to a particular fluorescent protein or translation inhibitor. Our results show that fluorescent-protein calibration can be condition-dependent and cannot be assumed to transfer across physiological perturbations.

DNA-based technologies rely on short, transient hybridization events, but selecting sequences with desired kinetic properties remains largely empirical because hybridization kinetics are difficult to predict from sequence and slow to measure one sequence at a time. Here, we introduce SPARXS-Hyb, an implementation of SPARXS (Single-molecule Parallel Analysis for Rapid eXploration of Sequence space) for multiplexed sequence-resolved screening of DNA hybridization. Using a surface-immobilized docking-strand library and a quencher-labelled imager-strand library, we screened 128 different DNA sequences in a single kinetic measurement, exposing all sequences to identical experimental conditions. This multiplexed approach removes a major confounding factor of serial measurements, allowing sequence-dependent differences to be compared directly. The resulting dataset reveals sequence-dependent transient binding behaviours and enabled us to identify a sequence with which an order-of-magnitude higher sampling rate can be achieved in DNA-PAINT (DNA points accumulation for imaging in nanoscale topography), a super-resolution microscopy technique based on DNA hybridization. By enabling multiplexed screening across a sequence library, SPARXS-Hyb provides a route to kinetics-guided sequence selection for programmable transient interactions in DNA nanotechnology.

Living cells are highly dynamic and densely crowded environments in which organelles such as vesicles undergo continuous motion that is essential for cellular processes. Therefore, accurate tracking of individual organelles is crucial for understanding intercellular dynamics and functions. However, precise tracking of individual organelles in living cells remains challenging due to high organelle densities, frequent particle overlap, and the coexistence of stationary and motile organelles. In particular, stationary organelles can obscure the trajectories of moving organelles, leading to tracking errors and fragmented tracks. To overcome these challenges, we developed Multiple Particle Tracking via Velocity Filtering (MPT-vVF), an unbiased, semi-automated tracking framework that incorporates a mathematically derived velocity-filtering algorithm to selectively identify and track moving organelles with high accuracy in crowded intracellular environments. MPT-vVF integrates denoising, background subtraction, and a velocity-matching detection step that discriminates true particle motion from noise based on spatiotemporal continuity, followed by robust trajectory linking. We demonstrate that MPT-vVF can accurately resolve nanometer-scale displacements of immobilized beads, highlighting its high tracking precision. We also validate the robustness of MPT-vVF by quantifying the transport of brain-derived neurotrophic factor (BDNF)-mRFP-containing vesicles in living hippocampal neurons. Furthermore, MPT-vVF reveals that exposure to 50-nm nanoplastics impairs vesicular transport, reducing both travel length and speed of BDNF-containing vesicles in living neurons. These findings establish MPT-vVF as a powerful method for quantitative analysis of intracellular organelles in crowded living cells and suggest its broad application to biophysics, cell biology, and soft matter research.

Transition-metal ion Forster resonance energy transfer (tmFRET) measures short-range distances and conformational changes in proteins, but its distance sensitivity depends on the donor-acceptor pair. Here, we establish the metal-chelating noncanonical amino acid bipyridylalanine (bpyA) as a flexible, genetically encoded acceptor platform whose Forster distance (R0) can be tuned without changing protein labeling. Cu2+, Cu2+ complexed with nitrilotriacetic acid (NTA), and Fe2+ form distinct acceptors at the same bpyA site, allowing distance sensitivity to be adjusted by changing the solution conditions. We incorporated bpyA at two sites in maltose binding protein (MBP) and labeled nearby engineered cysteines with the donor fluorophore BODIPY. The resulting [Cu(bpyA)]2+, [CuNTA(bpyA)]2+, and [Fe(bpyA)]2+ acceptors provided R0 values ranging from 13.9 to 21.2 [A]. Analysis of time-resolved tmFRET data recovered donor-acceptor distance distributions and maltose-dependent changes consistent with structural predictions. At the 237/295 site pair, both copper and iron acceptors resolved the maltose-dependent decrease in distance. At the 309/322 site pair, [Cu(bpyA)]2+ did not appreciably resolve the distance increase, whereas adding NTA shifted R0 toward shorter distances and enabled detection of an increase from 12.1 to 14.4 [A]. These results demonstrate that a single bpyA-labeled protein preparation supports multiple tmFRET acceptors, providing a cysteine-independent acceptor-labeling strategy with adjustable distance sensitivity.

Adhesion of cells is the key factor determining functioning of multicellular organisms. Viscoelastic properties of cells can be studied by multiple methods. However, attractiveness of cells or extracellular matrix without the elastic component (dispersive adhesion) is not accessible. We present an extension of force spectrometry technology: the Multivalent Adhesive Probe Atomic Force Microscopy (MAPA) that delivers dispersive adhesion maps of live cells and biosurfaces, and identifies differences unresolved by viscoelastic probing.

Cellular systems must act robustly to maintain organismal health, including maintaining biophysical properties that allow for appropriate cellular function, and adapting these properties through changes such as those that occur during cell division. However, we still lack tools to measure many of these physical properties with precision in the living cell. For example, the mechanical properties of the nucleoplasm, the fluid-like substance that fills the nucleus, have not been fully characterized. To investigate these properties, we have turned to the fission yeast Schizosaccharomyces pombe (S. pombe), a well-established, genetically tractable model organism that has been used extensively for studying a variety of cell biophysical processes and structures, including the cytoskeleton and cell division. It is an apt system for studying how the nucleus adapts over the course of the cell cycle, since it undergoes closed mitosis, where the nuclear envelope remains intact during cell division. Studying nucleoplasm properties over the course of closed mitosis may help reveal how nuclear volume, shape, surface area expansion, and chromosome segregation are linked and coordinated. To measure nucleoplasm material properties in S. pombe, we have paired Fluorescence Correlation Spectroscopy (FCS) with a photoswitchable fluorophore, enabling fine control over fluorescent intensity inside live cells. We infer material properties from FCS measurements, while the photoswitchable probe enables confocal imaging in conjunction with these measurements, yielding corresponding information about cellular state and dynamics. Interestingly, we find that nucleoplasm material properties do not vary significantly over the cell cycle. Future studies will use this tool to examine how diverse molecular and genetic perturbations alter nucleoplasmic properties, providing insight into how these properties maintain nuclear function and protect genomic integrity over the cell cycle and during development.

Fourier light field microscopy (FLFM) enables high-speed volumetric imaging by encoding multiple angular perspectives of a three-dimensional sample onto a single image. For this reason, FLFM is well-suited to sparse and rapidly evolving biological systems. To aid in the adoption of FLFM, we present OpenFLR, an open-source software framework for real-time volumetric reconstruction, three-dimensional particle tracking, and calcium image processing using FLFM. OpenFLR reconstruction is distributed as four interchangeable interfaces: a Python library, a command-line script, an interactive web application, and an ImageJ/micromanager plugin, so that the pipeline is accessible to both developers and bench biologists. Building on established Richardson-Lucy deconvolution, we use a hybrid experimental-computational PSF calibration strategy and a triangulation approach to tracking to extract particle positions in 3D directly from raw light field frames, bypassing reconstruction. We validate the complete pipeline on GCaMP6s recordings of freely behaving Hydra vulgaris, tracking sparse populations of neurons as they undergo large three-dimensional displacements.

Mechanical stretchers are used to physically expand biological samples for microscopic studies of mechanobiology. Existing stretchers suffer from limited strain capacity and non-quantifiable forces, along with focus drift and feature-tracking limitations during microscopy. We develop an iris-based cell stretcher that applies isotropic equibiaxial force to achieve aerial strain of 1664% that corresponds to 4.2x linear expansion for Tensile Expansion Microscopy (TExM), a super-resolution method which increases sample size above the diffraction limit of light. The TExM stretcher uses 3D-printed, printed circuit board (PCB) and cost-effective parts, is portable and automated. We verify the performance of the TExM cell stretcher using image analysis, achieving ~90% mechanical precision, a maximum of four degrees of arm angular deviation, precise speed control down to 0.01 cm/sec, and an average expansion resolution of 3.98 x 10^-3 x. Integrated strain gauge sensors confirm equal application of force by each arm of the stretcher throughout expansion and finite-element simulation and planar-polariscope photoelastic imaging characterize the uniform substrate stress distribution while indicating high stress at the substrate gripping area. An AutoTracking software in communication with the stretcher hardware and microscope enables continuous autofocus and feature tracking during TExM and fiducial markers verify uniform equibiaxial stretch. The stretcher is demonstrated with fixed NIH 3T3 fibroblasts and live HeLa cells, observing ~4x cellular expansion of fixed cell size and separation of live cell clusters, highlighting the potential of TExM cell stretcher for biological imaging.

Cardiac blood flow is a regulator of several important developmental and remodelling processes in the heart, including through fluid shear forces sensed by the endothelial cells lining the heart. However, optically mapping these flow fields in the complex 3D geometry of the heart is challenging even in transparent animal models such as the zebrafish. One of the main challenges is the difficulty in measuring the out-of-plane (axial) velocity component, preventing accurate mapping of the complete 3-component-3-dimension (3C-3D) blood flow velocity field; image-based techniques such as microscopic particle image velocimetry ({micro}PIV) traditionally only provide the in-plane flow components. Here we present a computational approach to achieve full time-varying 3C-3D blood flow vector mapping using a standard selective plane illumination microscope (SPIM), based on robust cardiac phase assignment, precise measurement-driven registration of sequentially acquired z-stacks, and PIV data fusion from multiple sample orientations. Our approach holds the key to understanding the complex dynamic flow fields within the developing heart, and their role in shaping cardiac development.

The high energy-efficiency of life is hard to understand only with classical physics. Many efforts have been made to study its mechanism based on quantum mechanics; the progress is nevertheless slow due to lack of experimental evidence with living cells. Here, combining experiments on cells, tissues and mitochondria with a theoretical model, we demonstrate a quantum state of mitochondria, which can be employed to modulate ATP production in living cells. We found an anomalous 71.0-THz oscillation mode only in living cells and tissues, which is highly determined by intact structure of mitochondria, and cannot be assigned to any specific molecules. Based on experimental data, a quantum model of light-matter coupling was introduced to trace the origin of this mode. Our calculations suggest a quantum superposition state of functional mitochondrion that forms by the coupling of light and lipid CH2 bonds in functional cristae, and induces a splitting of the intrinsic CH2 vibration mode of 87 THz to two levels at 71 THz and 103 THz, respectively. The former can be observed only in living cells and tissues; whereas the latter falls in the range (90-110 THz) of biomolecular and water vibrations, thus indistinguishable. Additional experiments revealed this mitochondrial quantum state able to serve as an efficient channel to modulate ATP production. Our findings provide a quantum mechanics view for understanding living cells, and it will be interesting to further explore whether such quantum state could act as a channel for energy metabolism, and even information transmission in life.